Dynamical Heterogeneity and the interplay between activated and mode coupling 

dynamics in supercooled liquids 
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We present a theoretical analysis of the dynamic structure factor (DSF) of a liquid at and below the 
mode coupling critical temperature Tc, by developing a self-consistent theoretical treatment which 
includes the contributions both from continuous diffusion, described using general two coupling 
parameter (F12) mode coupling theory (MCT), and from the activated hopping, described using the 
random first order transition (RFOT) theory, incorporating the effect of dynamical heterogeneity. 
The theory is valid over the whole temperature plane and shows correct limiting MCT like behavior 
above Tc and goes over to the RFOT theory near the glass transition temperature, Tg. Between Tc 
and Tg, the theory predicts that neither the continuous diffusion, described by pure mode coupling 
theory, nor the hopping motion alone suffices but both contribute to the dynamics while interacting 
with each other. We show that the interplay between the two contributions conspires to modify the 
relaxation behavior of the DSF from what would be predicted by a theory with a complete static 
Gaussian barrier distribution in a manner that may be described as a facilitation effect. Close to Tc, 
coupling between the short time part of MCT dynamics and hopping reduces the stretching given 
by the F12-MCT theory significantly and accelerates structural relaxation. As the temperature is 
progressively lowered below Tc, the equations yield a crossover from MCT dominated regime to 
the hopping dominated regime. In the combined theory the dynamical heterogeneity is modified 
because the low barrier components interact with the MCT dynamics to enhance the relaxation rate 
below Tc and reduces the stretching that would otherwise arise from an input static barrier height 
distribution. Many of these results can be explained from an analytical treatment of the combined 
equation of motion. 



I. INTRODUCTION 



In an earlier article we showed how to connect self- 
consistently the mode coupling theory (MCT) with the 
random first order transition theory (RFOT) to describe 
the dynamics of a liquid above and below the mode cou- 
pling transition temperature, Tc The resulting dy- 
namics includes both the diffusive dynamics described by 
MCT and the hopping dynamics described by RFOT the- 
ory. Although other earlier attempts to include both hop- 
ping and mode coupling dynamics within one theoretical 
scheme have been made [3, 01 : the merit of our calcula- 
tion was the use of hopping dynamics, determined via the 
RFOT theory, thus acknowledging in accord with exper- 
iments that the hopping rate decreases with the config- 
urational entropy. Because of the feedback between the 
structural relaxation (which includes contribution from 
both continuous and hopping dynamics) and the viscos- 
ity, hopping has a non-linear effect on the total dynam- 
ics. Due to the self-consistent nature of the calculation 
there was hopping induced softening of the growth of the 
frequency dependent viscosity with decreasing tempera- 
ture and this in turn helped the relaxation of the MCT 
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contribution to the structural relaxation. Thus, the the- 
ory predicts that below along with the input hopping 
dynamics there is an additional hopping induced continu- 
ous diffusion which was absent when hopping was frozen. 
The time scale of relaxation is thus found to be faster 
than that predicted by hopping motion alone because it 
now also includes the contribution from the continuous 
dynamics. This effect is key to showing explicitly that no 
strict localization transition takes place at Tc, in accord 
with long standing arguments [1, Q . 

To keep the theory analytically tractable in our earlier 
work we used a simpler version of the MCT which in- 
cluded only one coupling parameter A and neglected the 
distribution of barrier heights in the hopping dynamics 
predicted by RFOT theory As a result the a relax- 
ation was nearly exponential. In this present article in or- 
der to address the origin of the stretching of the long time 
a relaxation dynamics we examine not only a two cou- 
pling parameter MCT (Gotze's F12 model) having both 
Ai and A2 term which in combination can directly result 
in stretching [3| but we also incorporate the static barrier 
height distribution, from RFOT theory, in the hopping 
dynamics. The term containing Ai describes the coupling 
of the density relaxation to a static field (which may de- 
scribe a localized defect or a static inhomogeneity in the 
density of the system) and that containing A2 describes 
the self coupling. As shown by Gotze and co-workers, the 
F12 model which formally describes static inhomogeneity 
present in the system predicts a stretching of the a re- 
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laxation dynamics above Tc It is not clear precisely 
how such a static inhomogeneity would in fact be gener- 
ated above the microscopic Tc^ however, such a scenario 
is perhaps viable at temperatures below Tc but there the 
hopping dynamics must also contribute significantly. 

Computer simulation studies of atomic displacements 
in supercooled binary mixture systems strongly suggest 
the coexistence of continuous diffusion and hopping as 
mechanisms of mass transport @. These studies show 
that hopping events are often followed by enhanced con- 
tinuous diffusion. These two mechanisms can obviously, 
therefore, interact cooperatively with each other. 

The present analysis provides a quantitative descrip- 
tion of the non-linear interaction between continuous dif- 
fusion and dynamically heterogeneous activated hopping. 
It is shown that at and just below Tc, hopping helps un- 
lock the continuous diffusion which now becomes more 
effective than the hopping would be by itself. We further 
find that below Tc, the stretching of the relaxation com- 
bines the effects of activation inhomogeneity and static 
inhomogeneity. The barrier height distribution takes care 
of the dynamic inhomogeneity in the system and becomes 
the primary source of the stretching of the dynamics 
much below Tc, but the MCT effects play a role in the 
low barrier components to enhance the rate of short time 
diffusion. 

The organization of the rest of the paper is as follows. 
In the next section we describe the theoretical scheme. In 
section III we present several analytical results that can 
be derived for the combined theory. Section IV contains 
numerical results. Section V concludes with a discussion 
on the results. 



II. THEORETICAL SCHEME 

In our earlier article we showed that activated dynam- 
ics or hopping opens up an extra channel for the struc- 
tural relaxation The continuous dynamics was cal- 
culated using the one coupling parameter (A) MCT. In 
describing the activated motions the probability of a sin- 
gle hop was calculated from RFOT theory which con- 
nects the height of the free energy barrier to the con- 
figurational entropy. For simplicity we had considered a 
single value of the barrier height for each temperature 
although RFOT theory shows that this barrier is in fact 
distributed. The present theory uses the same scheme 
of calculation with some modifications to understand the 
relation to previous MCT efforts to cope with nonexpo- 
nential relaxation. (1) The one coupling parameter (A) 
MCT is extended to incorporate the two coupling pa- 
rameters (Ai, A2) MCT (the F12 model) which Gotze has 
used to address the effect of static inhomogeneity on the 
dynamics above Tc- (2) The single valued barrier height 
for the activated dynamics is replaced by a distribution 
of barrier heights in accord with RFOT theory. 

The previous article used two different mathematical 
schemes to combine the hopping and the continuous dy- 



namics (described by MCT)[l|. In one of the schemes 
the full intermediate scattering function was written as 
a product of a hopping and a MCT part using the sep- 
aration of timescales between the MCT dynamics and 
the hopping dynamics. In the second scheme the strict 
parallelism of hopping and continuous motion was more 
transparent. The structure of the equation in the sec- 
ond scheme is similar to that obtained by Gotze and 
coworkers Q and Das and Mazenko from more de- 
tailed microscopic derivations. Both the schemes give 
nearly identical results. This further adds credence to 
the first scheme. In the present paper, therefore, we 
will work with the first scheme (easier to implement and 
also in this scheme the continuous diffusion and hopping 
dynamics can be investigated separately) although the 
extensions made here can also be incorporated into the 
second scheme in a similar manner. 

The total intermediate scattering function can be writ- 
ten as, 



(j){q,t) ~ (l)MCT{q,t)(t>hopiq,t). 



(1) 



Here <t>MCT{q,t) is the MCT part of the intermediate 
scattering function, which is now self consistently calcu- 
lated with (j){q, t), and its equation of motion is given by, 

+ Xlill f dt'<i3{t')^MCT{t-t') 



(2) 



+ \2^l / dt'(f'{t')^MCT{t-t')^0 

Jo 

In the above equation the fourth term on the left hand 
side describes the coupling of (f>MCTiq,t) with a static 
field which is meant to describe the defects or the in- 
homogeneity in the system, according to Gotze 0. The 
fifth term on the left describes the coupling of 4>MCT{q, t) 
with itself (the self coupling term). Unlike the earlier 
model the present model contains two order parameters. 
In the absence of hopping, the MCT transition would 
now take place not at a single point but at many points 
on the Ai — A2 plane. 

In eqUthe hopping part of the intermediate scattering 
function is give by 4)hop{q,t)- The contribution from a 
single hopping event to the scattering function was de- 
rived in our earlier paper It can be written as, 



(l>lop{q) 



s + Khop{q) 



where, 
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FIG. 1: The (j){t) calculated from the unified theory, and the 
results of the idealized MCT are plotted for three tempera- 
tures. The solid lines A, B, C correspond to the full (j}{t) for 
e — —0.1, 6 = and e — 0.5, respectively. The dashed lines 
are the idealized MCT results. In this plot A — 0.568 which 
means that the MCT dynamics without hopping above Tc is 
stretched with stretching parameter Pmct = 0.8. The acti- 
vated dynamics is calculated with barrier height distribution 
such that /^f^-^*'" = 0.5. 



In the above expression of the hopping kernel, P is 
the average hopping rate which is a function of the 
free energy barrier height, AF and is given hy P — 
■^exp{—AF/kBT) The free energy barrier is calcu- 
lated from RFOT theory [§]. = f^i'^'^ is the region par- 
ticipating in hopping where ^ is calculated from RFOT 
theory. Vp is the volume of a single particle in the sys- 
tem. cIl is the Lindemann length. In this model kernel a 
typical hopping event involves an uncorrelated displace- 
ment of particles by a Lindemann length. More complex 
kernels that encode correlations between movements are 
also possible. 

Now if we consider a distribution of barrier heights 
then the contribution from multiple hoppings to the in- 
termediate structure function can be written as, 



(l>hop{q,t) 



bl,p{t)V{AF)dAF 



V{AF) is considered to be Gaussian. With a Gaussian 
distribution of barrier heights the relaxation function is 
known to fit well to a stretched exponential, where the 
stretching depends on the width of the Gaussian f?| . 

In figure 1 we plot the results of the unified theory and 
also that of the ideahzed MCT (F12 model). The plots 
are given for e = —.1, e = and e = 0.5 which correspond 
to temperatures above, at and below Tc respectively. The 
calculations are done at A = 0.568 which for idealized 
MCT above (for e < ) predicts a stretched relax- 
ation with the MCT stretching parameter, Pmct = 0.8. 
For the activated dynamics we have considered Gaussian 
distribution of barriers that would predict the stretching 



parameter for the hopping dynamics, P^op^'^ = 0.5. As 
expected, hopping does not have much effect above Tc 
but below Tc the unified theory continues to show struc- 
tural relaxation where the idealized MCT would have 
predicted strict localization transition. The longtime dy- 
namics is stretched at all the temperatures. But as will 
be discussed later the stretching parameter for the total 
dynamics is different from Pmct and P^^p^'^ ■ value 
depends on the interaction between the hopping and the 
MCT dynamics and changes with temperature. 



III. THE EFFECT OF HOPPING ON THE MCT 
DYNAMICS ABOVE, AT AND BELOW Tc 
:ANALYTICAL RESULTS 

In the earlier article we showed numerically that hop- 
ping has very little effect on the MCT dynamics above 
Tc and has a nonlinear effect on the MCT relaxation 
timescale at Tc. In this article we analytically investi- 
gate the effect of hopping on the MCT dynamics above, 
at and below Tc. For this study we first examine the 
one parameter MCT dynamics, i.e. we consider Ai = 0. 
This is a special case of the two parameter model. The 
MCT transition takes place at A2 = 4 for Ai = 0. Ini- 
tially we also take the barrier height distribution to be a 
delta function corresponding to a single hopping barrier. 
With these simplifications we use eqs[l][5]for the following 
analysis of the effect of hopping on the MCT timescale. 

When Ai = 0, Eql5]can be rewritten as. 



^l<j}MCT{z) 
1 + Z(f)MCT{z) 



= z 



-i^^LT\-,VLl\<^^(t)\{z) (6) 



where LT stands for Laplace transform. Since we are 
interested in the a relaxation timescale we will take the 
longtime limit of the above equation. In the longtime 
limit where z ^ 0, eqE] reduces to, 



— — — «7 + LT\2nc[(j) {t)\{z) 

1 -I- Z(j)MCT[Z) 



(7) 



In the above equation the second term on the right 
hand side is the self coupling term which makes dominant 
contribution at high density. But at low liquid density 
the self coupling term can be neglected and the solution 
of eqlT] is given by. 



9MCT 



(t) = e- 



-Kat 



(8) 



where Kq = is the inverse timescale of longtime decay 
in the normal liquid regime. 

Eq[7]in the time plane can be written as, 

1 ■ /■* 

(/>MCt(0 = -^(pMCrit) + A2 / dt'<j)^{t')(j)MCTit - <')(9) 
J<o Jo 

From the earlier analysis [11, Q we know that the solution 
of eq|5] is exponential in the longtime. So we simplify 
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4>MCT{t) using its longtime form, (j>MCT{t) = ae"^"'^'^*, 
where a is the prefactor. Presently we also take the hop- 
ping part of the intermediate scattering function to be 
an exponential, 4>hop{t) — e~'^''°p*. Using these expres- 
sions for 4>MCT and 4>hop in eq[9]we get an expression for 
Kmct in terms of Khop, Ko, X2 and a. 



Kmct — 2 



{2Khop-{l-a^X2)Ko 



(10) 



+ ^ {2Khop - (1 - a^X2)Koy + SKhopKo 
In the absence of hopping the above expression reduces 



to, 



K 



MCT 



= (1 - a^X2)Ko 



(11) 



Thus we see that the self coupling term leads to an in- 
creases in the timescale of relaxation when compared 
with the bare relaxation timescale K^^. EqfTT] further 
predicts that Kmct goes to zero or the relaxation time 
approaches infinity as 0^X2 approaches one. Thus the 
analysis shows that in the absence of hopping strict lo- 
calization takes place at a^A2 — 1. From previous studies 
we know at Tc, a = 1/2 and A2 = 4 [l|, [^, [l^], thus at 
T = TcT a?X2 indeed becomes unity. 

We will now analyze Kmct in the presence of hopping 
in three different regions, above Tc, at and around Tc and 
below Tr. 



A. Above Tc 

Above the mode coupling transition temperature Tc, 
it was shown by explicit calculation that the timescale 
of hopping dynamics is so much longer than direct re- 
laxation in the system that it can be neglected [l|. The 
expression of Kmct then reduces to. 



Kmct ^ {I - a'X2)K, 



(12) 



Thus in accord with our earlier numerical calculation [l| 
above the mode coupling transition temperature the hop- 
ping does not have significant effect on the a relaxation 
timescale. 



B. At and around Tc 

From the earlier studies we know that at the transition 
temperature, a^A2 = 1 [1, 9]. Thus cqlTT] reduces exactly 
to. 



Kmct = -Khop + ^ Kop + ^KhopKo (13) 

In the above expression if we take some reasonable value 
for Ko and Khop, we find that. 



Kmct — \/ 2KhopKo 



(14) 



Thus we find that Kmct has a nonlinear dependence on 
Khop which is also in accord with our earlier numerical 
calculations [1]. The coupling of the hopping dynamics 
with the short time part of the MCT dynamics (liquid 
like dynamics) leads to this nonlinear dependence. Thus 
we find that coupling between the short time part of the 
MCT dynamics and hoppings leads to an MCT part of 
the structural relaxation timescale which is much faster 
than the hopping timescale. This is a critical effect and is 
found at and near T^. This result corroborates our ear- 
lier findings that a single hopping event leads to many 
continuous diffusion events which are the primary means 
of structural relaxation in this region. In the next sub- 
section we will find that this scenario changes as we go 
lower and lower in temperature. 



C. Below Tc 

Much below the transition temperature, a^A2 >> 1. 
Eg I 111 can be rewritten as. 



Kmct — 



iia^X2 ~ \)Ko + 2Khop) 



,r \ SKhopK o 

■i'^^f^iia^X2-l)Ko 



2Khop) 



(15) 
'}] 



Since at low temperatures ((a^A2 — l)ifo + "^KhopY » 
^KhopKo thus we can write. 



Kmct 



^KhopKo 



{{0^X2 - l)Ko + 2Khop) 



(16) 



Now if we further take into consideration that (a^A2 
l)Ko » 2Khop then the above equation reduces to. 



Kmct 



2Kh 



op 



(a2A2 - 1) 



(17) 



Thus much below the transition temperature the 
timescale of the MCT dynamics and also the total dy- 
namics becomes slaved to the hopping timescale. Analyz- 
ing ea llTl an important observation can be made about 
the relaxation timescale. It is known that as we lower the 
temperature both A2 and a increases. Thus the denomi- 
nator in eq ll7l will increase as we lower the temperature. 
From the analysis in the earlier subsection we know that 
initially to start with, below Tc Kmct >> Khop, then as 
we keep lowering the temperature then depending on the 
temperature dependence of A2 and a ,Kmct — Khop- 
But as we further lower the temperature then slowly 
Kmct << Khop- Thus in this regime although there 
will be hopping induced continuous diffusion but the pri- 
mary mode of the structural relaxation becomes direct 
activated hopping itself. 
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IV. INTERPLAY BETWEEN HOPPING AND 
MCT DYNAMICS AT AND BELOW T = Tc : 
NUMERICAL RESULTS 

Continuing from our analysis where both MCT and 
hopping dynamics are assumed to be exponential, here 
we will present some numerical results for the general 
case where both hopping and MCT dynamics can be 
stretched. For these general cases we will try to under- 
stand the effect of hopping on the MCT dynamics. 

For the calculation we need to solve eql3] numerically. 
It is well known that due to the self consistent nature 
of the equation, its numerical solution becomes a nearly 
Herculean task around the mode coupling transition tem- 
perature. In the presence of hopping due to the disparate 
timescales present in the system and the convolution in 
eql3] which involves all these timescales, the time of cal- 
culation increases many fold. However the scheme pro- 
posed by Fuchs et al fll| allows the calculation to be done 
much faster. Both (j){t) and 4>MCT{t) vary more slowly 
for longer times than they do at short times. The essen- 
tial idea involved in the scheme is to separate the slow 
and the fast variables and treat them differently in the 
convolution. The short time part of (j)(t) and (j) Merit) 
are calculated exactly with very small stepsize and they 
are then used as input to carryout the calculation for the 
long time part of the same. We have exactly followed 
the scheme presented in reference [ll| except for a minor 
modification (one extra term) in eq.29 of reference [Tlj 
when the integration timestep is an odd number. In our 
calculations the total timestep N=1000. 

For this study we have picked three points on the Ai, 
A2 plane. In all the cases to understand the temperature 
effect the values of Ai, A2 were varied according to the 
expressions given bellow. 



_ (2A-1) 



A 



(1 + (1-A)2) 



A2 = + e 



A(l-A) 



A2 (1 + (1-A)2) 



(18) 



(19) 



where the value of A determines the value of Ai and A2 
ai T = Tc. e is a measure of the distance from the MCT 
transition temperature. As in our earlier model calcu- 
lation, 1\ the values of Q.* and 7* are kept unity and 
the scaling time as Ips. In the numerical calculations, 
to understand the temperature effect although we have 
changed the values of Ai and A2 but we have kept all the 
other parameters constant including the hopping barri- 
ers. For the calculation of the hopping part we have used 
eqlUwith a Gaussian distribution of barriers. According 
to the RFOT theory, the mean barrier height and thus 
the hopping timescale should change with temperature 
0, Q . In its simplest form without taking barrier soft- 
ening into consideration the mean barrier height can be 
written as, AF/keT = rnks/sc - ^'^^ 



Sc is the configurational entropy, Acp is the jump in the 
specific heat and Tk is the Kauzmann temperature 0, [1] . 
However in this present calculation to clearly assign any 
change in dynamics due to the change in Ai and A2 value 
we have kept the mean of the distribution fixed at about 
8.8 fcsT. As mentioned earlier the stretching in the hop- 
ping dynamics is determined by the width of the Gaus- 
sian distribution of the barriers. The broader the dis- 
tribution the more stretched is the dynamics. We have 
varied the width of the distribution such that the /3^*J^**'^ 
varies from 0.2 to 0.8. Along with the Pf^^p^'^ value the 
timescale of hopping also changes, which has been taken 
into account in our calculation. 

For all the cases the Ptotai values are plotted against 
Phop^'^ for different e values, where Ptotai is the stretching 
parameter for (j){t) and (3^^^'^'^ is the same for 4>hop{t). 



A. Case 1: A = 0.5 

In the first case we consider an example when A ~ 0.5. 
This value of A implies that at T = Tc, Ai = and A2 = 4 
and the MCT dynamics just above Tc is exponential. For 
this case we vary e from — 1. In figure 2, the Ptotai val- 
ues are plotted against (3'^^^'^'^ for the three e values. As 
expected for e = the j3totai — 1- This is because al- 
though without hopping structural relaxation is arrested 
but once hopping is present the relaxation is dominated 
by the MCT dynamics. As we increase the epsilon value 
we find that the effect of hopping is stronger and the dy- 
namics gets stretched. But we also notice that for very 
low P^^p^'^ value the dynamics is less stretched than the 
static barrier distribution would indicate. This is because 
lower Z?^*"**'^ value implies a broader barrier height distri- 
bution which means we populate both lower and higher 
barrier heights. The small barrier hoppings actually cou- 
ple to the liquid like part of the MCT dynamics (whose 
timescale is given by K'^^) and the relaxation is faster 
and dominated by this liquid like MCT dynamics which 
has a much shorter timescale. As we increase the e value 
MCT gets more slaved to hopping and thus the time scale 
of relaxation increases and the liquid like MCT dynamics 
becomes less important. If we further lower the temper- 
ature (or increase the e value) we would find that even 
for small /3^op**'^ values the total dynamics is determined 
primarily by hopping. 

The interplay of hopping with MCT nonlinearities can 
be thought of as a quantitative formulation of "facilita- 
tion effects" [7|. As Xia and Wolynes pointed out hop- 
ping events interact if they occur near each other. This 
is accounted for by the MCT nonlinearity. In the Xia- 
Wolynes treatment the corresponding effect led to the 
cutoff of the relaxation time distribution, on the slow 
side, owing to the renewal of a mosaic cell's environment 
through hops. This resulted in an increased (3 from that 
obtained from the static Gaussian model, as occurs here 
too [2]. 
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FIG. 2: The stretching parameter for the total structural re- 
laxation, 4>{t), Ptotai is plotted against the same for cf)hop{t), 
Phop^'^ for three different e values. The solid circle is for e = 0, 
the solid square is for e = 0.5 and the solid diamond is for 
e — 1. In this plot the A = 0.5 which means that the MCT 
dynamics without hopping above Tc is exponential. 



also find that for high /J^op**"^ values, as e increases, the 
dynamics gets more and more dominated by the hopping 
dynamics. Nevertheless, for small P^^p^'^ the scenario is 
a little different. For smaller e values the low barrier hop- 
pings couple to the liquid like part of the MCT dynamics 
and the dynamics is less stretched and also much faster. 
However, as we increase e, MCT dynamics gets more and 
more slaved to hopping dynamics and the MCT relax- 
ation timescale becomes proportional to the hopping re- 
laxation timescale (similar to that shown in Eq|17p . Thus 
the effect of the coupling of low barrier hoppings with the 
liquid like part of the MCT dynamics, on the total MCT 
dynamics, reduces. If we further lower the temperature 
then we would find that even for jS^^^"^"^ = 0.2 the total 
dynamics follows the hopping dynamics. 



The results in case 1 and case 2 look quite similar 
but more detailed observation reveals that the value of 
Phop^'^^ where Ptotai begins to increase, is smaller for case 
2 (where MCT dynamics itself is more stretched) than 
it is in case 1. As discussed earlier, the reason Ptotai in- 
creases for small P^op^'^ is that the small barrier hopping 
gets coupled to the liquid like part of the MCT dynam- 
ics allowing the total structure to relax. Now in case 2, 
the MCT dynamics is itself stretched thus the effect of 
the small barrier hopping on the MCT dynamics will be 
much less effective when compared to case 1. This trend 
becomes clearer when we study the next case where the 
MCT dynamics is much more stretched. 



Case 3: A = 0.75 



0.4 I 1 1 1 1 1 1 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 

FIG. 3: The stretching parameter for the total structural re- 
laxation, (j){t), Ptotai is plotted against the same for (f>hop{t), 
Phop"'" for three different e values. The solid circle is for e = 0, 
the solid square is for e = 0.5 and the solid diamond is for 
e = 1. In this plot the A = 0.568 which means that the 
MCT dynamics without hopping above Tc is stretched with 
stretching parameter, Pmct ~ 0.8. 



B. Case 2: A = 0.568 

In this case we consider that A = 0.568. This value of A 
implies that at T = T^, \i = .4215 and A2 = 3.09958 and 
the MCT dynamics, (without hopping) just above Tc, 
would already be stretched with (3mct = 0.8. For this 
case we vary e from 0—1. In figure 3, the Ptotai values are 
plotted against /^^o"*'^ for the three e values. The results 
are similar to those obtained for the first case. For e = 
the Ptotai — 0-8 which implies that the total dynamics 
is still determined primarily by the MCT dynamics. We 



In this case we consider that A = 0.75. This value of A 
implies that at T = T^ Ai = .889 and A2 = 1.778 and the 
MCT dynamics (without hopping) just above Tc would 
already be stretched with Pmct = 0.5. For this case we 
vary e from — 1. In figure 4, the Ptotai values are plot- 
ted against Pf^^p '''^ for the three e values. The results are 
similar to that obtained for case 1 and case 2 for e = 0. 
But for higher e values unlike in case 1 and 2, Ptotai con- 
tinuously decreases with Pf^^p^'^- This is because as dis- 
cussed before, since MCT dynamics is already stretched, 
the structural relaxation due to low barrier hopping is 
less effective. Also note that the Ptotai decreases with 
Phop^'^ value but it is neither equal to /S^^p**'^, nor equal 
to the pre-transition Pmct value. At these temperatures 
although MCT dynamics is slaved to hopping but both 
the channels of relaxation are almost equally effective. If 
we compare the Ptotai values for e = 0.5 and 1, we will 
find that in most of the cases Ptotai for e = 1 has a lower 
value. This is because at lower temperatures MCT dy- 
namics becomes a less effective relaxation channel. At 
further lower temperatures (higher e values) Ptotai will 
follow Pf^^^^'^ more closely. 
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FIG. 4: The stretching parameter for the total structural re- 
laxation, 4>{t), Ptotai is plotted against the same for cf)hop{t), 
Phop^'^ for three different e values. The solid circle is for e = 0, 
the solid square is for e = 0.5 and the solid diamond is for 
e = 1. In this plot the A = 0.75 which means that the MCT 
dynamics without hopping above Tc is stretched with stretch- 
ing parameter, Pmct = 0.5. 



V. CONCLUDING REMARKS 

MCT and the RFOT theory provide a unified the- 
ory of relaxation over the whole temperature plane 
Even without dynamical heterogeneity the theory suc- 
cessfully predicts the decay of the structural relaxation 
below mode coupling transition temperature, Tc, con- 
firming there is no strict localization transition at T^. 
Without dynamical heterogeneity of the instantons the 
coupled theory lead to an exponential a relaxation, but 
in the laboratory the a relaxation is generally stretched. 
In the present article we examined both a two parame- 
ter MCT model 0] and more realistically one that also 
included a barrier height distribution that gives rise to 
stretching in the hopping dynamics by itself [71]. 

The study has been carried out for different stretching 



parameters of the MCT dynamics, that is by changing A 
values in the F12 model and also for different stretching 
parameters of the hopping dynamics obtained by chang- 
ing the width of the distribution of barrier heights. 

To summarize, the main conclusions of the present 
work is that the continuous dynamics, described here 
within the F12— MCT formalism , and the activated hop- 
ping dynamics, described here using RFOT theory, in- 
teract in a non-linear fashion to give rise to dynamical 
features which are distinct from both. MCT by itself of 
course cannot describe dynamics below its critical tem- 
perature. We find that hopping facilitates the continuous 
dynamics channel and in the process the effects of hop- 
ping on the relaxation decreases. Thus, one finds the 
stretching parameter arising solely from distribution of 
hopping barrier energies in RFOT is increased by the 
mode coupling terms. It is also found that when MCT 
dynamics is less stretched then the effect of hopping on 
the MCT dynamics is less pronounced and the hopping 
dominated regime moves to a lower temperature. On the 
other hand, for more stretched MCT dynamics, due to 
the larger overlap of MCT and hopping timescales, hop- 
ping begins to dominate at a higher temperature. 

We have already mentioned the need for using a more 
complex hopping kernel than that given by Eqs.3-5. In 
particular, one needs to include the effects of mode cou- 
pling softening on the barrier height distribution. This is 
a feed-back effect of unleashing the mode coupling relax- 
ation channels due to hopping, on the barrier height dis- 
tribution itself. This non-linear feed-back is expected to 
shift the distribution to lower barrier heights and in turn 
accelerate mode coupling relaxation which can further 
enhance hopping. The whole system of equations needs 
to be solved self-consistently. To achieve this, we need to 
understand more quantitatively the effects of softening 
on the barrier height distribution. 
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